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SHOCK  TUBE  MEASUREMENTS  OE  ISO-OCTANE  IGNITION  TIMES 
AND  OH  CONCENTRATION  TIME  HISTORIES 

D.  F.  DAVIDSON,  M.  A.  OEHLSCHLAEGER, }.  T.  HEREON  and  R.  K.  HANSON 
High  Temperature  Gasdynamics  Laboratory 
Mechanical  Engineering  Department 
Stanford  University,  Stanford,  CA  94305,  USA 


Ignition  times  and  OH  radical  concentration  time  histories  were  measured  behind  reflected  shockwaves 
in  iso-octane/02/Ar  mixtures.  Initial  reflected  shock  conditions  were  in  the  ranges  1177  to  2009  K  and 
1.18  to  8.17  atm,  with  fuel  concentrations  of  100  ppm  to  1%  and  equivalence  ratios  from  0.25  to  2.  Ignition 
times  were  measured  using  endwall  emission  of  CH  and  sidewall  pressure.  OH  concentrations  were  mea¬ 
sured  using  narrow-linewidth  ring-dye  laser  absorption  of  the  Ri(5)  line  of  the  OH  A-X  (0,0)  band  at 
306.5  nm.  The  ignition  time  data  and  OH  concentration  time-history  measurements  were  compared  with 
model  predictions  of  four  current  iso-octane  oxidation  mechanisms,  and  the  implications  of  these  com¬ 
parisons  are  discussed.  To  our  knowledge,  these  data  provide  the  first  extensive  measurements  of  low-fuel- 
concentration  ignition  times  and  OH  concentration  time  histories  for  iso-octane  autoignition,  and  hence 
provide  a  critical  contribution  to  the  database  needed  for  validation  of  a  detailed  mechanism  for  this 
primary  reference  fuel. 


Introduction 

Iso-octane  is  an  important  component  of  practical 
fuels.  It  is  also  used  as  the  high-oetane-number  pri¬ 
mary  referenee  fuel  in  the  development  of  surrogate 
fuel  mixtures  for  gasoline  engine  ehemistry  model¬ 
ing.  Increasing  the  fraction  of  iso-octane  relative  to 
n-heptane  in  primary  referenee  fuel  mixtures  is  well 
known  to  reduee  knock  in  engines.  This  effect  is  at¬ 
tributed  to  the  strong  radieal  scavenging,  primarily 
of  OH,  of  the  iso-octane  oxidation  pathways  in  com¬ 
petition  with  the  reactive  n-heptane  oxidation  path¬ 
ways  [1].  Few,  if  any,  studies  exist  that  provide  direct 
information  on  the  OH  coneentration  time  histories 
that  occur  during  iso-octane  autoignition. 

The  chemistry  of  iso-oetane  oxidation  is  also  evi¬ 
dent  in  shock  tube  ignition  time  studies  where  the 
branehed-ehain  iso-octane  has  been  shown  to  have 
longer  ignition  times  than  the  straight-chain  n-oc- 
tane  [2].  However,  the  use  of  iso-octane  ignition 
times  as  model  validation  targets  is  eomplieated  by 
the  faet  that  there  does  not  appear  to  be  a  eonsensus 
among  the  various  shook  tube  studies  [2-6]  and  that 
few  studies  have  eorrelated  their  findings  with  those 
of  other  workers  to  allow  quantitative  comparison. 

The  chemical  kinetic  validation  of  iso-octane  re¬ 
action  mechanisms  currently  relies  on  a  wide  variety 
of  target  data  including  ignition-delay  times,  fiame 
speeds,  and  species  concentration  time  histories  in 
various  combustors  and  fiow  reactors.  Currently, 
only  stable  speeies  sueh  as  O2,  CO,  CO2,  and  H2O, 
and  long-lived  hydroearbon  intermediate  species, 
mostly  measured  using  gas  ehromatography,  are 


used  as  targets.  Very  limited  information  is  available 
on  the  important  transient  radieal  pool  (O,  H,  OH, 
HO2,  and  CH3)  responsible  for  H-abstraetion  reac¬ 
tions  or  branching  processes.  This  is  not  an  oversight 
on  the  part  of  these  mechanism  development  studies 
(see,  for  example.  Refs.  [7-10]);  these  data  simply 
do  not  exist  for  larger  fuel  molecules,  yet  they  are 
urgently  needed  if  the  predictive  ability  of  these 
meehanisms  is  to  be  tied  to  the  aetual  kinetics  of  the 
radical  pool  population. 

Many  of  the  eurrent  iso-oetane  reaction  mecha¬ 
nisms  are  derived  from  the  large  referenee  fuel 
meehanism  developed  at  Lawrenee  Livermore  Na¬ 
tional  Laboratory  (LLNL)  by  Westbrook  and  co¬ 
workers.  In  this  paper,  we  will  examine  three  of 
these  mechanisms:  Curran  et  al.  [7],  Pitsch  et  al.  [8], 
and  Davis  and  Law  [9],  as  well  as  a  mechanism  de¬ 
veloped  separately  in  Europe  by  Ranzi  et  al.  [10]. 

The  very  large  Curran  et  al.  iso-oetane  meehanism 
has  approximately  1000  species  and  4000  reactions 
and  represents  a  relatively  reeent  version  of  the  full 
LLNL  mechanism.  The  Davis  and  Law  mechanism 
(69  species  and  406  reactions)  was  formed  by  adding 
iso-oetane-specific  reactions  from  Curran  et  al.  to 
the  Held  et  al.  [11]  reduced  meehanism  for  n-hep- 
tane.  In  this  mechanism,  the  major  olefin  products 
of  iso-octane  decomposition  are  C3H6  and  i-C4H8, 
through  which  all  major  reaction  products  pass.  The 
Pitseh  et  al.  mechanism  includes  47  speeies  and  134 
reactions  and  is  an  extension  of  their  reduced  mech¬ 
anism  and  was  developed  for  use  in  full  engine  eal- 
culations  [8].  The  Ranzi  et  al.  mechanism  uses  145 
species  and  2500  reaetions  and  was  developed  using 
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Eig.  1.  Iso-octane  ignition:  CH  emission  profile  and  PZT 
pressure  trace.  Initial  shock  conditions:  0.25%  iso-octane, 
6.25%  O2  1269  K,  8.07  atm.  Ignition  time  is  952  jus. 

reaction  classes  and  a  lumped  kinetic  scheme  and 
includes  submechanisms  from  their  earlier  work  on 
n-heptane.  (See  Ref.  [10]  for  further  details.)  The 
need  to  form  smaller  mechanisms  by  reducing  the 
number  of  species  and  reactions,  and  therefore  mak¬ 
ing  these  mechanisms  more  amenable  to  use  in  full 
engine  calculations,  must  be  tempered  by  the  need 
to  accurately  capture  the  important  engineering 
characteristics  of  the  fuel.  The  ability  of  these  mech¬ 
anisms  to  do  this  in  the  case  of  high-temperature 
ignition  times  and  OH  concentration  time  histories 
will  be  demonstrated  in  this  study 

In  the  hydrocarbon  fuel  research  program  at  Stan¬ 
ford,  we  are  attempting  to  improve  and  extend  the 
database  of  kinetic  targets  used  in  the  testing  and 
validation  of  detailed  (and  reduced)  mechanisms. 
Recent  work  has  included  the  measurement  of  ig¬ 
nition  times  and  OH  and  C2H4  concentration  time 
histories  for  four  n-alkanes:  propane,  n-butane,  n- 
heptane,  and  n-decane  [12,13].  In  these  studies,  sig¬ 
nificant  differences  were  noted  between  the  mea¬ 
sured  OH  and  C2H4  concentration  time  profiles  and 
those  predicted  by  several  current  mechanisms.  The 
value  of  this  type  of  data  was  evident,  as  previously, 
there  was  no  other  direct  method  to  test  the  ability 
of  these  mechanisms  to  predict  radical  pool  popu¬ 
lations. 

In  this  paper,  we  present  new  data  for  ignition 
times  and  OH  concentration  time  histories  in  the 
iso-octane  oxidation  system.  These  data  are  corre¬ 
lated  and  compared  with  modeling  results  using  the 
four  kinetic  mechanisms  noted  above.  The  kinetic 
implications  of  these  comparisons  are  also  discussed. 

Experimental  Method 

All  ignition  times  were  measured  in  the  refiected 
shock  region  of  either  a  14.0  or  15.24  cm  diameter, 
high-purity,  helium-driven,  shock  tube.  The  shock 


velocity  at  the  endwall  was  calculated  from  four  in¬ 
cident  shock  wave  velocity  measurements  obtained 
over  the  last  1.5  m  of  the  shock  tube.  Refiected 
shock  conditions  were  determined  from  the  stan¬ 
dard  one-dimensional  shock  relations  and  the  Sandia 
thermodynamic  database,  which  was  amended  to  in¬ 
clude  the  thermodynamic  data  for  iso-octane  from 
Rurcat  and  McBride  [14].  Uncertainty  in  the  initial 
refiected  shock  temperature  is  estimated  at  ±  10  K, 
resulting  in  about  a  10%  uncertainty  in  the  measured 
ignition  time.  Research  grade  argon  and  oxygen  and 
99.7+%  HPLC  grade  iso-octane  (Aldrich)  were 
used  in  all  mixtures.  Initial  refiected  shock  condi¬ 
tions  of  these  measurements  were  in  the  ranges  of 
1177-2009  K  and  1.18-8.17  atm,  with  fuel  concen¬ 
trations  of  100  ppm  to  1%  and  equivalence  ratios, 
0,  from  0.25  to  2. 

Ignition  time  measurements  are  subject  to  a  high 
degree  of  experimental  scatter  and  uncertainty  if  the 
composition  of  the  test  mixture  is  not  accurately 
measured,  or  if  these  measurements  are  made  at  a 
significant  distance  from  the  endwall.  (See  Horning 
et  al.  [15]  for  a  discussion  of  these  effects  and  efforts 
made  to  minimize  such  error  in  our  facility)  Ignition 
times  in  this  study  were  measured  using  CH  emis¬ 
sion  viewed  through  the  shock  tube  endwall.  CH 
emission  was  collected  using  a  431  nm  filter  with  a 
spectral  width  of  10  nm  and  a  fast  (1.6  //s  risetime) 
silicon  photodetector.  The  endwall  ignition  time  is 
defined  as  the  time  interval  between  the  shock  ar¬ 
rival  at  the  endwall,  which  is  determined  from  the 
incident  shock  velocity,  and  the  maximum  rate  of 
increase  in  the  CH  emission  signal  (see  Fig.  1). 

OH  absorption  data  were  acquired  across  the  di¬ 
ameter  of  the  shock  tube  2  cm  from  the  endwall 
using  narrow-linewidth  ring-dye  laser  absorption. 
The  absolute  OH  concentrations  were  measured  us¬ 
ing  a  standard  two-beam  laser  absorption  technique 
developed  in  our  laboratory  (see  Herbon  et  al.  [16] 
for  a  complete  description)  using  the  well-charac¬ 
terized  Ri(5)  line  of  the  OH  A^i:+  -  X^n{0,0) 
band.  The  largest  uncertainties  in  the  OH  concen¬ 
tration  stem  from  the  OH  oscillator  strength  ( ±  2%) 
and  the  collisional  broadening  parameter,  yielding 
total  uncertainties  in  Xqh  in  the  range  of  2%-5%  for 
the  conditions  of  the  current  work.  Off-line  mea¬ 
surements  revealed  evidence  of  a  weak,  but  broadly 
absorbing  species.  This  absorption,  likely  from  con¬ 
jugated  olefins,  is  similar  to  that  seen  in  our  earlier 
study  of  OH  absorption  in  JP-10  ignition  [17],  but 
not  seen  in  our  study  of  n-heptane  and  other  n-al- 
kanes  [12].  The  interference  absorption  was  sub¬ 
tracted  from  the  on-line  absorption  data  to  deter¬ 
mine  the  corrected  OH  concentration  time  histories. 
An  example  OH  data  trace  and  the  results  of  the 
subtraction  process  are  shown  in  Fig.  2.  The  large 
signal-to-noise  ratio  and  high  sensitivity  of  this  di¬ 
agnostic  permit  measurement  of  parts  per  million 
levels  of  OH  at  microsecond  time  resolution. 
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Eig.  2.  OH  absorption:  on-line,  off-line,  and  eorreeted 
profile.  Initial  shoek  eonditions:  1.0%  iso-oetane,  6.25%  O2, 
1656  K,  1.34  atm.  Refleeted  shoek  arrival,  t  —  0;  ignition, 
t  =  116  JUS. 


Eig.  3.  Variation  of  ignition  time  with  temperature: 
0  —  1.0.  Initial  shoek  eonditions:  0.5%  iso-oetane,  6.25% 
O2  1.3  atm.  Eilled  squares,  eurrent  study;  thin  solid  line, 
fit  to  data  with  “  43.9  keal/mol;  solid  line,  Davis  and 
Law  [9];  dashed  line,  Ranzi  et  al.  [10];  dot-dashed  line, 
Pitseh  et  al.  [8];  dot-dot-dashed  line,  Curran  et  al.  [7]. 

Experimental  Results  and  Discussion 


Eig.  4.  Variation  of  ignition  time  with  fuel  eoneentration. 
Initial  shoek  eonditions:  0  —  1.0,  1500  K,  1.3  atm.  Eilled 
squares  and  thin  solid  line,  eurrent  study;  eross-square, 
Niemitz  et  al.  [4];  open  triangle,  Vermeer  et  al.  [3];  open 
eirele,  Bureat  et  al.  [2];  solid  line  (all  data  normalized  to 
1.3  atm  by  Davis  and  Law  [9];  dashed  line,  Ranzi 

et  al.  [10];  dot-dashed  line,  Pitseh  et  al.  [8];  dot-dot-dashed 
line,  Curran  et  al.  [7]. 


Eig.  5.  Variation  of  ignition  time  with  equivalenee  ratio 
for  fixed  O2  eoneentration.  Initial  shoek  eonditions:  6.25% 
O2  1500  K,  1.3  atm.  Eilled  squares  and  thin  solid  line,  eur¬ 
rent  study;  solid  line,  Davis  and  Law  [9];  dashed  line,  Ranzi 
et  al.  [10];  dot-dashed  line,  Pitseh  et  al.  [8]. 


Ignition  Time  Data 


The  ignition  time  data  are  presented  as  follows. 
First,  the  effects  of  temperature  and  fuel  concentra¬ 
tion  on  ignition  time  are  presented,  for  a  fixed  equiv¬ 
alence  ratio  of  1  (Figs.  3  and  4).  Second,  the  effects 
of  equivalence  ratio  on  ignition  time  are  shown  (Fig. 
5).  Third,  a  correlation  is  presented  for  the  ignition 
time  data  set  (Fig.  6). 

The  current  ignition  time  data,  in  general,  have 
very  small  scatter,  as  a  result  of  accurate  mixture  con¬ 
centration  and  shock-speed  determinations,  and  the 
use  of  endwall  CH  emission  as  a  direct  indicator  of 
ignition  time.  In  the  correlation  presented  below, 


there  is  strong  evidence  that  the  ignition  times  scale 
as  and  this  normalization  is  applied  to  all  the 

data  of  Figs.  3-6. 

The  1.3  atm,  0.5%  data  (shown  in  Fig.  3)  have  an 
activation  energy  of  43.9  keal/mol,  in  relatively  good 
overall  agreement  with  all  four  models.  All  four 
models  predict  ignition  times  within  a  factor  of  two 
of  the  experiment.  It  should  be  noted  that  the  cur¬ 
rent  calculations  were  performed  using  a  constant 
volume  constraint  and  nitrogen  as  the  carrier  gas  (as 
these  mechanisms  were  optimized  for  air.)  Small  var¬ 
iations  can  exist  between  model  predictions  with 
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Eig.  6.  Correlated  ignition  times:  various  shock  tube 
studies.  All  data  normalized  to  1.3  atm  and  1.68%  iso-oc- 
tane,  21%  O2  using  exp(  —  232Xfuei).  Eilled  squares, 

current  study;  cross-squares,  Niemitz  et  al.  [4];  open  tri¬ 
angles,  Vermeer  et  al.  [3];  open  circles,  Burcat  et  al.  [2]; 
open  diamonds,  Nixon  et  al.  [5]. 

M  =  Ar  and  M  =  N2  because  of  the  differences  in 
heat  capacity  of  Ar  and  N2  (but  this  effect  is  likely 
to  be  small  because  of  the  low  concentrations  of  re¬ 
actants  used  in  these  experiments)  and  because  of 
differences  in  the  collision  efficiencies  of  Ar  and  N2 
for  dissociation  reactions.  However,  the  effects  of 
these  differences  are  minor  on  our  calculations  for 
Tign-  (See  Ref.  [18]  for  a  study  confirming  this  ob¬ 
servation  in  CH4  oxidation.) 

The  variation  of  ignition  time  with  fuel  concentra¬ 
tion  presented  in  Fig.  4  deviates  from  a  simple 
power-law  dependence.  This  is  particularly  evident 
at  the  highest  concentrations  where  the  present 
study  shows  a  strong  decrease  in  ignition  time  as 
concentration  increases.  This  was  also  seen  in  the 
Burcat  et  al.  [2]  and  the  Niemitz  et  al.  [4]  data,  but 
was  not  seen  in  the  Vermeer  et  al.  [3]  data.  At  low 
to  intermediate  concentrations,  the  Ranzi  et  al  and 
the  Davis  and  Law  models  closely  follow  the  data, 
while  the  Pitsch  et  al.  and  Curran  et  al.  models  pre¬ 
dict  shorter  ignition  times  by  approximately  a  factor 
of  two.  At  high  concentrations,  the  four  model  pre¬ 
dictions  fall  within  the  scatter  of  the  three  previous 
studies  [2-4]  and  the  current  data. 

The  variation  of  the  ignition  time  with  equivalence 
ratio  (presented  in  Fig.  5)  shows  a  relatively  weak 
effect  on  ignition  time  of  fuel  concentration  for  rich 
mixtures  and  fixed  oxygen  concentration.  At  1500  K 
and  1.3  atm,  as  the  equivalence  ratio  varies  from  I.O 
to  2.0,  the  measured  ignition  time  only  varies  29% 
from  400  to  515  jus.  The  Davis  and  Law,  Pitsch  et 
al.,  and  Ranzi  et  al.  mechanisms  all  predict  similar 
variation  of  the  ignition  time  with  equivalence  ratio 
for  lean  mixtures,  but  do  not  predict  the  rolloff  for 
rich  mixtures  seen  in  the  experiments. 

The  variation  of  ignition  time  with  pressure  and 
equivalence  ratio  can  be  scaled  with  some  accuracy 


assuming  a  simple  power  law  [15].  However,  the 
rapid  decrease  in  ignition  times  with  increased  fuel 
concentration  suggests  that  the  variation  of  ignition 
time  with  fuel  concentrations  could  be  better  cor¬ 
related  with  an  exponential  function. 

Using  these  relationships,  a  correlation  for  our  95 
ignition  time  data  points  was  found: 

Tign  =  4.50  X  exp(-232Xf^ei)^^-®^ 

X  exp(44,780/K[eal/mol/K]r)  (1) 

with  ignition  time  T^gn  in  seconds,  pressure  P  in  at¬ 
mospheres,  temperature  T  in  degrees  Kelvin  fuel 
mole  fraction  and  equivalence  ratio  0.  Data 
from  the  present  work  and  from  several  other  shock 
tube  studies  are  normalized  (for  differences  in  P  and 
Xfuei)  using  this  correlation  and  are  presented  in 
Fig.  6.  Here,  the  data  of  Niemitz  et  al.  are  most 
consistent  with  the  present  study,  and  the  ignition 
time  data  of  Burcat  et  al.  appear  to  converge  with 
the  current  study  at  higher  temperatures.  The  data 
of  Vermeer  et  al.  converges  with  the  current  data  at 
lower  temperatures,  but  diverge  markedly  at  higher 
temperatures. 

The  highest  quality  shock  tube  data  are  obtained 
at  very  low  fuel  and  oxidizer  concentrations,  and  at 
short  test  times.  With  longer  test  times  (>1  ms), 
shock  tube  data  should  be  corrected  for  temperature 
and  pressure  changes  resulting  from  increased 
boundary  layer  growth.  With  higher  fuel  concentra¬ 
tions,  the  temperatures  and  pressures  that  occur 
during  the  rapid  release  of  energy  during  ignition 
are  not  well  characterized.  The  current  low-concen¬ 
tration  data  do  not  require  significant  correction,  and 
should  provide  an  accurate  determination  of  ignition 
time  over  the  conditions  studied. 

The  accurate,  low-scatter  data  in  Figs.  3-5  and  the 
correlation  presented  in  equation  I  should  be  useful 
in  testing  the  ability  of  current  iso-octane  mecha¬ 
nisms  to  predict  changes  in  ignition  time  with  pres¬ 
sure,  fuel  concentration,  equivalence  ratio,  and  tem¬ 
perature.  For  further  evaluation  of  these  detailed 
models,  data  on  radical  species  concentration  time 
histories  are  needed.  These  data  is  presented  in  the 
next  section. 

OH  Laser  Absorption  Data 

Representative  OH  concentration  time  histories 
(corrected)  are  shown  in  Fig.  7,  and  a  comparison  of 
one  time  history  with  modeling  from  several  mech¬ 
anisms  is  shown  in  Fig.  8. 

Table  I  summarizes  specific  details  of  the  OH  con¬ 
centration  time-history  profiles  (corrected  for  inter¬ 
ference  absorption).  In  this  table  (and  illustrated  in 
Fig.  8),  Tx(peak),  T2(min),  and  z^{50%)  are  the  times 
for  the  OH  mole  fraction  to  reach  the  first  peak 
value,  the  minimum  after  this  peak,  and  50%  of  the 
second  plateau  level,  respectively;  X(Ist  max), 
X(min),  and  X(2nd  max)  are  the  OH  mole  fraction 
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Eig.  7.  OH  absorption  data  (corrected).  Initial  condi¬ 
tions:  500  ppm,  iso-octane,  0.625%  O2.  Upper  trace, 
1736  K;  lower  trace,  1449  K. 

initial  peak,  the  minimum  or  plateau  value,  and  the 
postignition  peak  value,  respectively,  is  the  ab¬ 
sorption  coefficient  that  was  used  to  convert  the 
measured  absorbances  to  OH  mole  fraction  [16]. 
13(50%)  can  be  related  to  the  ignition  time.  Model 
calculations  indicate  that  13(50%)  is  nearly  coinci¬ 
dent  with  the  maximum  of  the  CH  concentration 
(which  is  closely  equivalent  to  the  maximum  rate  of 
change  of  the  CH  endwall  emission)  that  occurs  dur¬ 
ing  the  most  rapid  formation  of  radicals  in  the  final 


Eig.  8.  OH  absorption  data  (corrected)  and  comparison 
with  models.  Initial  shock  conditions:  500  ppm  iso-octane, 
0.625%  O2  1511  K,  1.51  atm.  Heavy  solid  line,  measure¬ 
ment;  solid  line,  Davis  and  Law  [9];  dashed  line,  Ranzi  et 
al.  [10];  dot-dashed  line,  Pitsch  et  al.  [8];  dot- dot- dashed 
line,  Curran  et  al.  [7]. 

stages  of  ignition.  The  measured  values  of  13(50%) 
are  in  good  agreement  with  the  simultaneously  mea¬ 
sured  endwall  ignition  times  if  the  small  sidewall 
time  correction  of  Petersen  et  al.  [18]  is  applied.  This 
correction  is  generally  required  for  measurements  of 
ignition  time  where  the  observation  is  made  away 
from  the  endwall. 


TABLE  1 

Summary  of  iso-octane  OH  absorption  data 


Ts 

P5 

K 

Ti  (peak)*  X  (1st  max) 

T2  (min)* 

X  (min) 

T3  (50%)* 

X  (2nd  max) 

(K) 

(atm) 

(atm“^  cm“^) 

(jus) 

(ppm) 

(jus) 

(ppm) 

(us) 

(ppm) 

0.05%  iso-octane,  0.625%  O2 

1736 

1.399 

132 

5 

44 

27 

36 

137 

456 

1614 

1.463 

151 

5 

30 

86 

9 

388 

352 

1511 

1.511 

173 

6 

22 

2 

1152 

288 

1449 

1.593 

185 

8 

16 

1 

2182 

245 

0.25%  iso-octane,  6.25%  O2 

1570 

1.412 

160 

159 

42 

62 

Sat.^ 

1467 

1.482 

178 

3 

93 

14 

12 

226 

Sat. 

1394 

1.533 

198 

4 

49 

49 

2 

621 

Sat. 

1299 

1.625 

209 

6 

16 

1 

1825 

Sat. 

0.5%  iso-octane,  6.25%  O2 

1507 

1.435 

154 

100 

7 

280 

Sat. 

1641 

1.373 

148 

149 

12 

34 

62 

Sat. 

1408 

1.495 

190 

3 

42 

2 

880 

Sat. 

1330 

1.554 

205 

4 

17 

4 

1870 

Sat. 

1.0%  iso-octane,  6.25%  O2 

1322 

1.576 

206 

3 

12 

1 

Sat. 

1398 

1.491 

192 

31 

1 

1422 

Sat. 

1506 

1.418 

172 

26 

3 

589 

Sat. 

1656 

1.336 

147 

103 

17 

12 

127 

Sat. 

^Sat.  =  transmission  ~  0. 
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CHEMICAL  KINETICS  IN  COMBUSTION— Kinetics  Measurements  and  Models 


As  can  be  seen  in  Fig.  7,  all  OH  absorption  data 
show  evidenee  of  the  same  coneentration  time  his¬ 
tory  structure.  At  early  times,  there  is  a  rapid  for¬ 
mation  of  OH  to  a  peak  X(lst  max)  simultaneous 
with  the  initial  decomposition  of  the  fuel  molecule. 
The  OH  coneentration  then  falls  to  a  well-defined 
plateau  X(min),  during  the  induction  period.  The 
scavenging  effeet  of  iso-oetane  and  its  deeomposi- 
tion  products  on  the  radical  population  is  clearly  evi¬ 
dent  in  the  suppressed  levels  of  OH  that  oeeur  in 
this  plateau  X(min).  This  suppression  is  especially 
evident  when  compared  to  similar  measurements 
done  in  n-heptane  [12]  where  the  OH  plateau  level 
before  ignition  does  not  drop  below  its  initial  for¬ 
mation  peak.  In  the  rapid  final  ignition  process,  OH 
eoneentration  rises  to  the  postignition  plateau  X(2nd 
max)  and  then  rounds  smoothly  to  a  seeond  plateau. 
This  second  plateau  level  falls  slowly  over  the  re¬ 
corded  millisecond  times  scale  of  these  experiments 
as  slower  three-body  recombination  processes  take 
place. 

The  predietions  of  the  four  mechanisms  are  com¬ 
pared  with  one  example  OH  coneentration  time  his¬ 
tory  data  trace  in  Fig.  8.  All  four  mechanisms  cap¬ 
ture  the  general  form  of  the  OH  coneentration 
profile.  The  Ranzi  et  al.  meehanism  most  closely 
models  the  strong  suppression  of  the  OH  radicals 
before  ignition,  but  both  the  Ranzi  et  al.  and  the 
Curran  et  al.  meehanisms  do  a  reasonable  job  at 
modeling  the  full  profile,  predicting  within  a  factor 
of  two,  the  initial  OH  peak.  The  smaller  Davis  and 
Law  and  Pitsch  et  al.  mechanisms  overpredict  the 
minimum  OH  by  faetors  of  3  and  10,  respeetively, 
and  the  first  peak  by  a  factor  of  3,  though  the  Davis 
and  Law  meehanism  still  does  an  excellent  job  at 
predicting  the  ignition  times. 

Sensitivity  analysis  using  the  Davis  and  Law  meeh¬ 
anism  for  the  conditions  of  Fig.  8  indicate  that  at 
early  times  the  formation  of  the  first  OH  peak  is 
controlled  by  the  reactions:  H  +  O2  ^  O  +  OH, 
iso-octane  ^  UC4Hg  +  i-C4Hg,  iso-oetane  ^  C7H45 
+  CH3,  and  i-C3H7  ^  H  +  C3H5;  the  subsequent 
OH  removal  is  eontrolled  by:  i-C3H7  ^  CH3  + 
C2H4.  As  expected  for  this  branched  alkane,  the  fate 
of  OH  radieals  is  intimately  related  to  the  ehoiee  of 
decomposition  pathways  and  rates  for  iso-octane, 
and  improvements  to  the  modeling  of  this  early  time 
OH  are  likely  possible  by  adjusting  the  rates  of  these 
pathways  within  their  uneertainty  limits. 

During  the  OH  plateau  minimum,  contribution 
analysis  indicates  that  there  is  a  competition  be¬ 
tween  OH  formation  by  H  +  O2  and  OH  removal 
by  reaetion  with  H2,  CH3,  CH2O,  C2H5,  and  other 
intermediate  product  species.  The  OH  formation 
rate  is  sensitive  to  the  fate  of  H-atoms  and  therefore 
shows  a  strong  sensitivity  to  reactions  such  as  H  + 
C3H5,  and  H  +  C3H5  that  remove  H-atoms,  and 
HCO  +  O2  that  lead  to  formation  of  H-atoms.  Im¬ 
provements  in  modeling  the  plateau  minimum  may 


come  from  the  inclusion  of  a  fuller  complement  of 
H-abstraetion  reaetions  for  the  intermediate  olefins 
and  OH  removal  reactions  for  all  intermediate  spe- 
eies,  as  well  as  adjustment  of  eritieal  rate  eoeffieients. 

While  ignition  time  and  OH  data  help  identify  key 
reaetions,  the  next  steps  in  meehanism  improvement 
are:  (I)  adjustment  of  critical  rate  eoeffieients,  within 
their  uneertainties;  (2)  direet  studies  of  key  reaetion 
rates  where  they  are  uneertain;  (3)  measured  time 
histories  of  other  speeies,  for  example,  CH3,  to  fur¬ 
ther  constrain  the  models;  and  (4)  comparative  in¬ 
vestigations  of  other  branehed  alkanes,  for  example, 
iso-butane  and  iso-pentane,  which  would  provide  in¬ 
formation  on  related  kinetic  structure. 

Conclusions 

A  new  database  of  ignition  times  and  OH  radical 
concentration  time  histories  was  generated  for  iso¬ 
octane  ignition.  A  comparison  of  the  ignition  time 
data  with  a  variety  of  eurrent  iso-oetane  oxidation 
mechanisms  revealed  that  these  mechanisms  repli¬ 
cate  some  of  the  features  of  the  ignition  time  data 
reasonably  well,  but  do  not  accurately  capture  the 
temperature  dependenee  of  the  ignition  time  or  the 
rolloff  behavior  with  rich  equivalence  ratios.  As  well, 
caution  is  advised  when  seleeting  modeling  targets 
ignition  times  at  high  coneentrations,  because  not  all 
the  published  data  are  consistent  with  the  present 
study.  A  comparison  of  the  measured  OH  time  his¬ 
tories  with  several  reeent  models  reveals  that  wide 
variations  occur  in  the  predictions  of  the  OH  con¬ 
eentration  time  histories  espeeially  for  the  smallest 
redueed  mechanisms.  Researchers  engaged  in 
model  development  should  find  these  new  data 
helpful  in  refining  their  models  and  when  used  in 
conjunetion  with  similar  data  for  n-heptane  per¬ 
formed  in  our  laboratory  [12]  should  provide  a  new 
database  for  primary  reference  fuel  meehanism  val¬ 
idation. 
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COMMENTS 


Frederick  L.  Dryer,  Princeton  University,  USA.  In  the 
eomparison  with  your  shoek  tube  ignition  delay  data  with 
those  of  other  investigators,  did  you  eonsider  that  tem¬ 
perature  eorreetions  of  the  ealeulated  temperature  might 
be  of  different  magnitudes  in  the  different  experiments? 

Do  you  have  any  meehanistie  or  experimental  explana¬ 
tion  for  the  observation  of  the  initial  OH  “spike”  found  in 
experiments  and  not  shown  by  any  of  the  meehanisms? 

Authors  Reply.  We  did  not  attempt  to  apply  any  tem¬ 
perature  eorreetions  to  our  data  or  to  the  published  data 
of  other  groups  when  we  eompared  them.  If  shoek  tubes 
of  substantially  smaller  diameters  are  used  for  ignition  time 
measurements,  there  may,  however,  be  some  reduetion  in 
the  apparent  aetivation  energy  due  to  boundary-layer- 
indueed  temperature  inereases  that  oeeur  in  shoek  tubes 
at  longer  test  times. 

Ror  a  very  short  duration  {S  jus)  during  the  passage  of 
the  shoek  front  aeross  the  laser  beam  path,  the  shoek  front 


momentarily  steers  the  diagnostie  beam  away  from  the  de- 
teetor  surfaee.  This  generates  a  sharp  oseillation  in  the  ab¬ 
sorption  data  at  time  zero.  This  spike  has  no  ehemieal  ki- 
netie  signifieanee. 

• 

William  J.  Pitz,  Lawrence  Livermore  Laboratory,  USA. 
In  our  work  on  the  autoignition  of  iso-oetane/air  mixtures 
under  eonditions  related  to  argone  knoek,  we  also  found 
that  the  OH  eoneentration  was  suppressed  prior  to  hot  ig¬ 
nition.  We  determined  through  ehemieal  kinetie  modeling 
that  the  formation  of  iso-butane  from  iso-oetane  and  the 
reaetion  of  iso-butene  with  OH  removes  OH  from  the  sys¬ 
tem.  This  reaetion  leads  to  resonantly  stabilized  iso-butenyl 
radieals  and  ehain  termination. 

Authors  Reply.  Thank  you  for  your  eomment. 


